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ABSTRACT 



We report the detection of strong and compact molecular line emission (in the 
CO J=3-2,4-3,6-5,7-6, 13 CO J=3-2, HCN and HCO+ J=4-3 transitions) from 
a cometary-shaped object (Carina-frEGGl) in the Carina star-forming region 
(SFR) previously classified as a photoevaporating protoplanetary disk (prop lyd). 
We derive a molecular mass of 0.35M o for Carina-frEGGl, which shows that 
it is not a proplyd, but belongs to a class of free-floating evaporating gas glob- 
ules (frEGGs) recently found in the Cygnus SFR by Sahai, Morris & Claussen 
(2012). Archival Adaptive Optics near-IR (Ks) images show a central hourglass- 
shaped nebula. The derived source luminosity (about 8 — 18 L Q ), the hourglass 
morphology, and the presence of collimated jets seen in HST images, imply the 
presence of a jet-driving, young, low-mass star deeply embedded in the dust in- 
side Carina-frEGGl. Our results suggest that the true nature of many or most 
such cometary-shaped objects seen in massive SFRs and previously labelled as 
proplyds has been misunderstood, and that these are really frEGGs. 

Subject headings: Stars: formation, Stars: pre-main sequence, Stars: protoplanetary 
disks, ISM: jets and outflows, ISM: individual objects: PCYC1173 
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1. Introduction 

Massive stars have a profound effect on stars that subsequently form in their vicinity, 
thus significantly influencing the initial mass function (IMF), the star formation efficiency, 
and the total extent and mass of star clusters. In massive star-forming regions (SFRs), the 
ionizing radiation from OB stars can produce photo-evaporating protoplanetary disks or 
proplyds, the most famous of which lie in the Orion Nebula (e.g., O'Dell, Wen, & Hu 1993; 
Ricci, Robberto, & Soderblom 2008). This radiation can also progressively ionize dense 
knots of nearby molecular material, forming evaporating gaseous globules, or EGGs (Hester 
et al. 1996). 

A new class of tadpole-shaped objects was recently discovered near the Cygnus OB2 
association (Sahai, Morris & Claussen 2012 [SMC12], Wright et al. 2012 [Wetall2]) - the 
most massive young association within 2 kpc of the Sun, located inside the Cygnus-X 
giant molecular cloud complex where considerable star formation is ongoing. The strong 
morphological similarity between these objects and proplyds led Wetall2 to conclude that 
these objects are "a unique class of photoevaporating, partially-embedded, young stellar 
objects" and unlikely to be EGGs, from a comparison of the large fraction (about 70 %) of 
these objects that appear to contain one or more young stars with the much lower value 
for M16 (15%). However, SMC12 concluded that Wetall2's argument is not compelling, as 
several unconsidered variables, including density and mass of parent cloud, and strength of 
radiation field and winds, can affect this fraction. In addition, the sizes of these objects 
(20,000-110,000 AU) are huge compared to the Orion proplyds (40-350 AU: Henney & 
O'Dell 1999). Recognizing the importance of these issues, SMC12 carried out molecular-line 
observations that reveal dense molecular cores associated with two of these objects (the 
Tadpole and the Goldfish), with total masses of cold (~ 10 — 15 K) molecular gas exceeding 
1-2 M Q , strongly favoring the EGG hypothesis. These data also conclusively rule out the 
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proplyd hypothesis since a proplyd is highly unlikely to harbor such a substantial mass of 
molecular gas. To date, no dense molecular medium has yet been found towards any known 
proplyd, and the low Orion proplyd disk masses (0.003-0.07) M & (Mann & Williams 2010), 
derived from continuum observations of dust emission, imply even smaller masses for the 
circumstellar material being evaporated from them. Finally, the evaporated material is 
likely to be primarily atomic, not molecular. 

SMC12 concluded that these tadpole-shaped objects are dense, star-forming molecular 
cores that originated in the Cygnus cloud and are now being photoevaporated by ultraviolet 
radiation from the CygOB2No. 8 cluster located ~10pc to their North- West, and shaped 
by ram pressure of strong, nearby wind sources. Continuing observations of 13 CO and 12 CO 
lines, as well as high-density tracers such as the J=3-2 lines of HCO + , HCN, HNC, and 
N 2 H + in a sample of about 20 such objects have strongly supported the SMC12 hypothesis 
(Sahai et al. 2012b). 

The discovery of these free-floating Evaporating Gas Globules (hereafter frEGGs) 
raises a fundamental question for star and protoplanetary disk formation in massive SFRs: 
what fraction of the previously classified proplyd objects are really EGGs? The distinction 
between frEGGs and proplyds is quite important, because in the former case, star formation 
and protostellar accretion are still very much under way, whereas proplyd disks are in much 
later stages in which accretion no longer has a strong effect on the protostellar evolution, 
although planet formation may still be in progress. We have therefore begun a study of 
the molecular gas content of objects previously classified as proplyds with the goal of 
determining their true nature. Smith, Bally & Morse (2003: SBM03) found many objects 
in the Carina SFR from ground-based optical emission-line images that closely resemble the 
proplyds found in Orion, and concluded that the former were proplyds as well in spite of 
their significantly larger (factor ~ 5) sizes. In this Letter, we report single-dish molecular 
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line observations of one such representative object from the set of 12 reported in SBM03, 
104632.9-600354 (hereafter Carina-frEGGl), and show that our data rule out the proplyd 
hypothesis for this object and support its frEGG nature. 



2. Observations & Results 

We observed Carina-frEGGl (Fig.QJ in the CO J=3-2, 4-3, 6-5 and 7-6 lines, the 13 CO 
J=3-2 and the HCN and HCO + J=4-3 lines, and the continuum at 350 and 870/im using 
the Atacama Pathfinder Experiment (Glisten et al. 2006) 12-m telescopqj during June 
23-26, 2012. Pointing was checked on nearby 77 Car, and found to be generally accurate 
within 3". 

The dual-color receiver FLASH+ (Heyminck et al. 2006) with wideband (4 GHz 
width) sideband-separating (2SB) SIS mixers was used to map the CO J=4-3 and 3-2 
lines simultaneously in a regular raster with 10" increments. Two XFFT spectrometers, 
each with a bandwidth of 2.5 GHz, but with 64 x 10 3 (32 x 10 3 ) channels for the low (high 
frequency) bands, were used to process a total IF band of 4 GHz (with 1 GHz of overlap). 
The weaker 13 CO(3-2), HCN and HCO + (4-3) lines were observed towards the nominal 
center position only. 

The dual-color heterodyne array receiver CHAMP+ (Kasemann et al. 2006), providing 
2x7 beams, was used to map the CO J=6-5 and 7-6 lines simultaneously, over a 40" x 40" 
region that was oversampled with 4" spacing (with all 7 beams of each sub-array covering a 
given grid position repeatedly). Two FFT spectrometers, each with a bandwidth of 1.5 GHz 
and 1024 channels, were connected to the individual beams, processing a total IF band of 



X APEX is a collaboration between the Max-Planck-Institut fuer Radioastronomie, the 
European Southern Observatory, and the Onsala Space Observatory 
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2.6 GHz (with 400 MHz overlap). The observations took place during excellent weather 
conditions with zenith precipitable water vapor, PWV, of 0.7 (0.4) mm in the low (high) 
frequency band. 

The spectra were taken with position-switching against an absolute reference position. 
Calibration was performed regularly every 10-15 min with a cold liquid nitrogen (LN2) 
load and an ambient temperature load. The data were processed with the APEX 
real-time calibration software (Muders et al. 2006). Beam sizes (FWHM: full- width at 
half-maxiumum) and main beam coupling efficiencies are given in TableHJ 

We obtained continuum observations at 350 (870) fim using the bolometer array 
receivers LABOCA (SABOCA) (Siringo et al. 2009, 2010), with a spatial resolution of 
7.8" (19") on 2012 Sep 10 (Sep 11) and PWV ~ 0.5 (0.3 - 0.4) mm. Data reduction was 
done with the BOA software (Schuller et al. 2009) following standard procedures, including 
iterative source modeling. 

The CO maps show the presence of a compact molecular globule (Fig.EJ), that appears 
unresolved at our highest angular resolution of 7.7" in the CO J=7-6 line. The line profiles 
(Fig.[3ti) show a central core and weak wings, extending to about ±6 km s" 1 from the line 
center as seen in our highest S/N profile (J=4-3), and suggesting the presence of an outflow. 
Gaussian fits to the line profiles result in a width (FWHM) of about 3.1km s" 1 in the 3-2 
and 4-3 lines, and 3.5 and 4.1km s _1 in the 6-5 and 7-6 lines, respectively, implying an 
increase in the wing contribution relative to the core for the higher excitation lines. The 
J=4-3 and 3-2 maps show the presence of low-level extended background or foreground 
emission around the compact source seen in the J=6-5 and 7-6 lines, but the contribution 
is small (<10%). 

In order to determine relative line intensity ratios, we convolved the CO 6-5 and 7-6 
maps to the 4-3 map resolution, 13.1" (we chose not to convolve these maps to the larger 
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3-2 beamsize because of possible contamination from the extended cloud to the East). 
Since the source is unresolved, the source's average brightness temperature in each line, 
Tr, depends on its emitting area. Conservatively assuming that all CO lines come from 
the dark region within the (ionized) bright periphery seen in the optical image, that has 
an area of 11.7arcsec 2 , we derive the radiation temperature, Tr, for each of the observed 
CO lines (Table 1). Using the online tool for the non-LTE radiative transfer code, RADEX 
(Van der Taak et al. 2007) to model these, we find that the 4-3 and 3-2 data can be fitted 
with a column density N(H 2 )~ 10 22 cm -2 and kinetic temperature ~ 40 K (assuming 
a CO/H2 abundance ratio fco ~ 10 -4 ), but the 6-5 and 7-6 line intensities predicted by 
this one-temperature model are weaker than observed (Fig.|3b), even with a model density 
high enough (rin 2 ~ 10 6 cm~ 3 ) to thermalize these lines, implying the presence of gas with 
Tfc >> 40 K in the globule. The "excess" emission in the 6-5 and 7-6 lines requires a hotter 
(> 150 K) optically-thin component assuming it has the same emitting the cooler 

component. The observed CO 6-5, 7-6 and HCO + 4-3 line intensities imply the presence 
of high gas densities, > 5 x 10 5 cm~ 3 . 

We derived the molecular mass in Carina- frEGGl from its 13 CO J=3-2 flux, using 
Eqn. (5) of Thi et al (2001), who used it to estimate protoplanetary disk gas masses 
from observations with the JCMT 15-m telescope (14" beam). We modified this equation 
to account for the APEX-12m's larger beam-size. We estimate a "beam- averaged" 13 CO 
J=3-2 optical depth, r 13 = 0.56 from the 12 CO/ 13 CO line intensity ratio, assuming 
T ex ( 13 CO)=T cx ( 12 CO)=40K. Assuming a 13 CO/H 2 abundance ratio f 13CO = 10" 6 and 
distance 2.3 kpc (e.g., Smith et al. 2004), we find a molecular mass of M mo \ = 0.35 M & . 

A 196 ± 65 mJy continuum source was detected at 350 /im, and an upper limit (40 mJy, 
la) obtained at 870 /im. Assuming the 350 /xm flux is due to optically-thin, thermal 
emission from dust at the same temperature as we derive above for the bulk of the gas (i.e., 
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Td = Tfc ~ 40 K), a gas-to-dust ratio of 100, and a dust opacity k(350 /mi) = 10cm 2 g -1 
(e.g., Miettinen et al. 2012), we derive a total mass of ~ 0.1 M Q . The discrepancy between 
this value and the mass derived from CO emission is not significant, considering the 
uncertainties in the adopted values of Td, k(350 /im) and fuco- 

3. Additional Multi- wavelength Data 

We have compiled multi-wavelength data on Carina-frEGGl using optical (HST), 
near-IR (2MASS and ESO) and mid-IR (Spitzer) data archives. Carina-frEGGl was 
imaged with ACS in Ha as part of a large-scale survey of the Carina SFR (Smith, Bally & 
Walborn 2010), and shows a dark tadpole- shaped globule (of length ~19,000 AU) with a 
bright periphery. The HST Ha (+[NII]) image of Carina-frEGGl shows a narrow, highly 
collimated jet (Fig.[Tb) - Smith et al. (2010) discuss this jet (HH1006) and infer that the 
putative central star related to the jet is not detected, presumably as a result of being 
deeply embedded within dust in the object. 

The 2MASS images show a faint red compact source in the J, H, and Ks bands. Using 
two nearby field stars to register the ACS and 2MASS Ks images, we find that the source 
is located at R.A. = 10 h 46 m 32 s .97, Dec = -60°03 m 53 s .5, roughly in the middle of the 
dark waist. We also retrieved archival Ks images towards Carina-frEGGl taken with the 
NACO+CONICA instrument on the ESO-VLT-U4 on 10 Jan 2007, with a total field-of-view 
27.8" x 27.8". Using standard IRAF packages, we generated a flat-field from these images 
using median-averaging and applied it to each of the 17 30-sec dithered exposures, which 
were then registered to a common reference frame and averaged. The angular resolution in 
the final image is about 0"1 as measured from radial intensity cuts of several neighbouring 
field stars. We find an extended faint hourglass-shaped nebula with a bright central region 
(Fig. [lb) with its symmetry axis aligned with that of the jet. The northern lobe is much 
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brighter than the southern one, implying that the hourglass structure is tilted such that its 
northern lobe is closer to us. In the central region of the northern lobe, which has a wide-U 
shape with its apex located at the center of the hourglass (inset, Fig.[Tb), the brightest 
feature is a small spur (of length ~50 AU) oriented at a position angle, PA ~ 12°. The spur 
is located slightly west of the jet axis and is not aligned with the latter (which is oriented 
at PA 10°). 

The Spitzer IRAC images (3.6-8 /im) show an elongated morphology similar to 
Carina-frEGGi's optical shape (but at much lower angular resolution). As in the Ks 
image, the shape is bipolar, with the northern half significantly brighter than the southern 
one. No central star can be seen in the Ks image or the IRAC images, implying that the 
former is deeply embedded within a dusty disk with very high line-of-sight extinction. In 
longer- wavelength Spitzer images (i.e., > 24/im), no isolated source can be detected at the 
location of Carina- frEGGl due to bright extended emission from foreground or background 
clouds in its vicinity. 

We have determined the SED of Carina-frEGGl from 1.25 to 870 /im. We carried 
out aperture photometry on the J and H-band 2MASS images (since the 2MASS catalog 
gives only a Ks magnitude). Aperture photometry for IRAC (channel 1 to 4) and MIPS 
(channel 1) is provided by Povich et al. (2011: source PCYC1173 in their catalog). Since 
no compact source can be seen in the MIPS 24 /im image at Carina-frEGGi's location, we 
treat Povich et al.'s measured 24 /j,m magnitude of 3 as an upper limit. The SED (Fig.HJ) 
continues to rise towards 350 /im - the longest wavelength for which we have a detection 
- indicating a substantial amount of emission from cool circumstellar dust in the far-IR 
range. For the source luminosity, we integrate the SED up to 350 fim to estimate L > 8 L Q ; 
a rough upper limit of 18 L & is provided by assuming flux values equal to the upper limits 
at 24 and 870 /xm. The observed fluxes have been corrected using an interstellar extinction 
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Ay = 1.5, that is intermediate between the value (1.78) given by Smith et al. (2004) and 
that (1.3) computed from a numerical algorithm provided by Hakkila et al. (1997), which 
estimates the Ay of a galactic source from its longitude, latitude and distance. 

We modeled the SED using an online tool that applies least-squares fitting to find 
the best models from a large set of pre-computed YSO models having accretion disks 
and infall envelopes with biconical outflow cavities (Robitaille et al. 2007). The ten 
best-fit models (Fig.HJ) have relatively low stellar effective temperatures (2660-3800 K), 
stellar masses (0.1-0.65 M Q ) and total luminosities (stellar + disk-accretion: 6-10 L Q ) and 
small ages (1100-3000 yr). These results support our inference of a deeply-embedded 
YSO in Carina-frEGGl, but must be treated with some caution. Although the inferred 
inclination of the outflow to the line-of-sight, i = 18°, together with the large circumstellar 
extinction along the disk's mid-plane (Ay > 100), is qualitatively consistent with the strong 
asymmetry in the brightness of the bipolar lobes in the Ks image, the projected opening 
angle of the cavity, 6 pro j = tan -1 (tan 9 C / sin i) = 23° (8 C ~ 7°, is the intrinsic opening angle 
in these models) appears smaller than observed. 



4. Discussion 



Our detection of a substantial mass of molecular gas in Carina-frEGGl shows that it 
is an frEGG, not a proplyd (as classified previously). Since the discovery of proplyds in 
Orion, many studies have reported the finding of proplyds in other SFRs with massive stars 
(e.g., Balog et al. 2006, Koenig et al. 2008, Brandner et al. 2000). However, in light of this 
study and that in SMC12, it is likely that the true nature of many or all of these objects 
has been misunderstood, and that some (or even all) of the previously classified proplyds in 
Carina, especially those which are significantly larger in size than the Orion proplyds, are 
really frEGGs. 
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Both proplyds and EGGs are unique probes of the effects of the harsh UV radiation 
and strong stellar winds from massive stars on the formation of lower-mass stars in their 
vicinity, thus it is of fundamental importance to be able to distinguish between these two 
classes. EGGs are most likely the surviving high density concentrations in a cloud as 
the ionization front sweeps through it. The formation of a protostar in the EGG may 
be induced by the compression of the latter by the high-pressure environment of the 
ambient HII region, supplemented by the even higher pressure at the high-density ionization 
front at the EGG surface. Photoevaporation sculpts out and exposes an EGG while the 
stellar object (s) in it is (are) still accreting mass, ultimately freezing the protostellar mass 
distribution at a relatively early stage in its evolution (e.g., Hester et al. 1996). FrEGGs 
like Carina-frEGGl (where the presence of collimated jets, central hourglass structure, and 
low luminosity, all imply the presence of a jet-driving, young, low-mass star) or the Tadpole 
(SMC12) represent the early formation phase of low or intermediate-mass stars within a 
massive SFR. Proplyds may represent the endpoint of evolution of such frEGGs, or may 
form coevally with the OB stars in these SFRs. Whether other objects like Carina-frEGGl 
are proplyds or EGGs has major implications for the star- format ion rate and the IMF in 
the massive SFRs where they are found. Outstanding questions that future research should 
address are: (i) do proplyds represent the evolutionary endpoint of frEGGS and (ii) is 
the formation of stars that may form in frEGGs induced by the surrounding HII region's 
compression or were such stars already in the process of formation in the relatively dense 
cloud cores that are susceptible to becoming frEGGS? 
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Table 1: Molecular Line Observations of Carina-frEGGl 
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a Main-beam line- intensities 

b Source brightness temperature derived by applying beam-dilution corrections 
c Derived from peak position in map convolved to 13.1" resolution 
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Fig. 1. — Images of Carina-frEGGl at different wavelengths (a) Spitzer IRAC 8/im, (b) 
HST F658N narrow-band (F658N filter: Ha+[NII]), and (c) VLT/NACO Ks (2//m). The 
panels are 14.6" x 24.4", with North up (as shown). Inset at right-middle shows a schematic 
of the hourglass- shaped central nebula (two possible curves are shown for the southern 
lobe's poorly-defined western periphery); inset at bottom- left corner shows an enlarged view 
(1.89" x 1.96") of the central source. 
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Fig. 2. — Map of Carina-frEGGl in CO J=6-5 emission integrated over the central 
5 km s _1 of the line profile, obtained with the APEX-12m. White circle shows the 8.7" beam- 
size (FWHM). 
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Fig. 3. — (a) Molecular line emission from Carina-frEGGl observed with the APEX-12m. 
The line profiles have been rescaled and shifted vertically for clarity. Weak absorption 
features in the CO J=4-3, 6-5 and 7-6 linse at Vi sr ~ —11 to —16 km s" 1 are due to 
emission in the reference position used for the position-switching. Intensities have not been 
corrected for main-beam efficiencies, (b) Observed radiation temperatures of CO lines from 
Carina-frEGGl (filled black squares) with 20% error bars, and a two-component model fit 
(open black squares). The lower temperature, optically-thick component at 40 K (red) fits 
only the 3-2 and 4-3 intensities; a second, hotter (~ 150 K), optically-thin component (blue) 
is needed to fit the low and high-J line intensities simultaneously. 
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Fig. 4. — Spectral energy distribution of Carina-frEGGl. Errors in the photometry are taken 
to be ±15% for data from 1.24 /zm to 8 /xm, and ±30% at 350 /im; black triangles show upper 
limits. The black curve shows the best-fit YSO model assuming ISM extinction Ay = 1.5 
and distance 2.3 kpc, and the gray lines show subsequent good fits. The dashed line shows 
the stellar photosphere from the best-fit model in the absence of circumstellar dust (but 
including ISM extinction). 



